home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / slaruv.z / slaruv
Encoding:
Text File  |  2002-10-03  |  3.1 KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAARRRRUUUUVVVV((((3333SSSS))))                                                          SSSSLLLLAAAARRRRUUUUVVVV((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLARUV - return a vector of n random real numbers from a uniform (0,1)
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE SLARUV( ISEED, N, X )
  13.  
  14.          INTEGER        N
  15.  
  16.          INTEGER        ISEED( 4 )
  17.  
  18.          REAL           X( N )
  19.  
  20. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  21.      These routines are part of the SCSL Scientific Library and can be loaded
  22.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  23.      directs the linker to use the multi-processor version of the library.
  24.  
  25.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  26.      4 bytes (32 bits). Another version of SCSL is available in which integers
  27.      are 8 bytes (64 bits).  This version allows the user access to larger
  28.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  29.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  30.      only one of the two versions; 4-byte integer and 8-byte integer library
  31.      calls cannot be mixed.
  32.  
  33. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  34.      SLARUV returns a vector of n random real numbers from a uniform (0,1)
  35.      distribution (n <= 128).
  36.  
  37.      This is an auxiliary routine called by SLARNV and CLARNV.
  38.  
  39.  
  40. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  41.      ISEED   (input/output) INTEGER array, dimension (4)
  42.              On entry, the seed of the random number generator; the array
  43.              elements must be between 0 and 4095, and ISEED(4) must be odd.
  44.              On exit, the seed is updated.
  45.  
  46.      N       (input) INTEGER
  47.              The number of random numbers to be generated. N <= 128.
  48.  
  49.      X       (output) REAL array, dimension (N)
  50.              The generated random numbers.
  51.  
  52. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  53.      This routine uses a multiplicative congruential method with modulus 2**48
  54.      and multiplier 33952834046453 (see G.S.Fishman,
  55.      2**b: an exhaustive analysis for b = 32 and a partial analysis for b =
  56.      48', Math. Comp. 189, pp 331-344, 1990).
  57.  
  58.      48-bit integers are stored in 4 integer array elements with 12 bits per
  59.      element. Hence the routine is portable across machines with integers of
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAARRRRUUUUVVVV((((3333SSSS))))                                                          SSSSLLLLAAAARRRRUUUUVVVV((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      32 bits or more.
  75.  
  76.  
  77. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  78.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  79.  
  80.      This man page is available only online.
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.